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The tribology of a sliding elastic continuum in contact with a disordered substrate is investigated 
analytically and numerically via a bead-spring model. The deterministic dynamics of this system 
exhibits a depinning transition at a finite driving force, with complex spatial-temporal dynamics 
including stick-slip events of all sizes. These behaviors can be understood completely by mapping 
the system to the well-known problem of a directed-path in higher- dimensional random media. 
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Pinning phenomena occur widely in nature and are 
important to many technological applications [1]. For 
instance, the pinning of flux lines suppresses dissipation 
in a Type-II superconductor, and the pinning of con- 
tact lines controls the spreading of a liquid on a solid. 
These are examples of driven elastic systems embed- 
ded in higher-dimensional random media. There have 
been much research activities in this area in recent years: 
The macroscopic properties of these systems, hereforth 
referred to as type (i) systems, have been elucidated 
through a combination of the functional renormalization- 
group (FRG) analysis [2,3] and numerous experimen- 
tal [4] and numerical [5-7] studies. In the extreme 
nonequilibrium limit where thermal fluctuations can be 
neglected and the dynamics is deterministic, one finds 
that a driving force exceeding a critical value is neces- 
sary to depin the system. A dynamic phase transition 
occurs at the depinning threshold, where the system ex- 
hibits complex stick-slip motion with avalanches of all 
sizes. The complex dynamics occurring in type (i) sys- 
tems results from the intricate balance of elasticity and 
random pinning forces near the onset of motion. Several 
universality classes classified by symmetries have been 
unveiled [8,9]. 

Another type of pinning phenomena which have at- 
tracted much interest is the friction which occurs between 
two sliding bodies interacting via a contact surface. This 
is exemplified by the Burridge-Knopoff model of earth- 
quake fault dynamics [10-12]. Other examples include 
boundary layer lubrication [13], and recent experiments 
on forced motion of elastic continuum (rubber and gel) 
over sticky substrates [14]. This type of systems, referred 
to here as type (ii) , also exhibit complex spatio-temporal 
dynamics, in the limit of very slow macroscopic motion. 
However, the dynamics of these systems are not as well 
understood. Recent numerical studies [11,12] of type (ii) 
systems relie on certain (chaotic) instabilities of the de- 
terministic dynamics to generate complexities. A sys- 
tematic description of spatio-temporal chaos is of course 
notoriously difficult. 

Since the interfaces of sliding bodies are almost always 
rough and irregular, and quenched randomness is known 



to generate complex dynamics (in type (i) systems), it is 
worthwhile to investigate the possible role of randomness 
in the dynamics of type (ii) systems. In this letter, we 
report analytical and numerical studies of a determinis- 
tically driven elastic system in contact with a disordered 
substrate. Our work is inspired by recent experimen- 
tal studies of friction of elastic continuum (stretched la- 
tex membranes and gel) on sticky substrates [14]. We 
model the inhomogeneous elastic medium by an array of 
beads connected by random springs. For simplicity, we 
focus on the properties of the ID version of such a ran- 
dom bead-spring array (i.e., a random chain), in contact 
with a disordered substrate (see Fig. 1). The effect of 
the 2-dimensional nature of the actual contact surface 
and the transverse dimension of the 3D elastic medium 
can be addressed by straightforwardly generalizing the 
results of the current study [15]. In our model, we assume 
that there is a large normal force confining the chain to 
the substrate, thus suppressing the transverse degrees of 
freedom [16]. We show below that our model of a ran- 
dom chain on a random substrate has a simple contin- 
uum limit, which we conjecture to be that of a directed 
path (DP) in 2-dimensional random media. The DP is 
the prototypical example of a type (i) systems for which 
much are known [1,17]. We demonstrate numerically that 
both the low temperature static properties of the random 
chain and its zero-temperature driven dynamics near the 
depinning threshold are indistinguishable from that of a 
randomly pinned DP. 
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FIG. 1. The model of a random chain being pushed across 
a disordered substrate. 

We study a one dimensional chain of beads coupled to 
their nearest neighbors by harmonic springs of random 
equilibrium lengths (see Fig. 1). The chain interacts with 
a (independently) random substrate described by the po- 
tential V(x), with V = 0, V(x)V(x') = A v S(x - x'). 
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Here, the overbar denotes disorder average. The (poten- 
tial) energy of the chain with N beads is given by 

N 

E [{r n }} = £ ( r «+i " r " " ««) 2 + V M , (1) 

n=l 

where r n denotes the position of the n th bead. The elas- 
ticity constant k is taken to be equal for all springs. 
The spring lengths a n are randomly chosen from a nar- 
row and bounded distribution with mean a, so that 
a n = a(l + Tj n ), Ti n <C 1, with Tfo = 0, Tiwq^ = A a 8 n , n ' ■ 

We consider deterministic and purely dissipative dy- 
namics of the system driven by a uniform external force 
F along the chain direction. In the strongly damped 
limit, the equations of motion is simply given by gradi- 
ent descent, i.e., 

xp 

H^d t r n {t) = -—+F, (2) 

where po is a kinetic coefficient. 

As we will soon show, the chain is completely pinned by 
the random potential at small driving forces. As the force 
increases above some threshold value, F c , the chain starts 
to move with a nonzero average velocity, v ~ (F — F C )P , 
for F > F c in the steady state. As is usual in critical 
phenomena, the behavior of the system near the (depin- 
ning) transition is characterized by a diverging correla- 
tion length, £ 3> a. It then becomes useful to coarse grain 
the discrete system (2). This will allow us to systemat- 
ically filter out the microscopic details and focus on a 
simpler continuum equation which preserves all features 
relevant to the critical phenomena. The coarse-grained 
density and currents are defined as 

r-x + a/2 

p(x,t)=a- 1 dx'Y^Hx' ~r n (t)), (3) 

Jx-a/2 n 
r-x + a/2 

j(x,t)=a~ 1 dx'J2 d trnS(x'-r n (t)). (4) 

Jx-a/2 n 

It is convenient to introduce a displacement-like field 
u(x,t), defined as p(x,t) = po (1 — d x u), where po = l/a 
is the average bead density. Then the continuity equation 
dtp(x,t) + d x j(x,t) = becomes d t u = j(x,t)/po. Using 
the equation of motion (2) in Eq. (4), then applying the 
coarse graining procedure and keeping the most relevant 
terms [15], we obtain: 

p^dtu = K d 2 x u - Pq 1 v8 x u + f[u, x] + F, (5) 

where K = na 2 is the chain stiffness, v(x,t) = 
Po(K d 2 u + f[u,x] + F) is the instantaneous velocity 
field as defined from j(x,t) = v(x,t)p(x,t), and 

/ = -d x V + 2p V P0 (x) sin [27rp (x - u - >y(x - u))] (6) 

results from the coarse-graining of the pinning force 
—dV/dr n - In (6), the sine term reflects the discrete 



nature of the bead-spring system [15], V Po denotes the 
Pg h -Fourier component of V, and 'y(s) = J s ds'ri(s') orig- 
inates from the random spring lengths. 

For a homogeneous chain (for which r) = 0), Eqs. (5) 
and (6) become similar to the well known random- 
amplitude sine-Gordon equation studied in the context of 
charge-density waves (CDW) [2], vortex array [18], and 
crystal growth [19]. The discrete translational symme- 
try u — > u + integer reflects in this case the invariance 
of the underlying discrete system upon a shift in the in- 
dex n labeling the beads. Compared to the usual CDW 
equation of motion, Eqs. (5) and (6) contain in addition 
the u- independent driving force d x V and the convective 
term linear in d x u. These terms have recently been in- 
troduced on phenomenological grounds [20,21], and their 
effects have been investigated in the limit of strong drive. 

The introduction of random springs leads to a factor 
7(x — u) which modulates the period of the sine term and 
hence breaking the discrete translational symmetry. This 
results from the non-uniformity of springs which destroys 
the relabeling invariance of the system [22]. The form of 
the probability distribution for f[u,x] is obviously com- 
plex. However, its leading moments are straightforward 
to compute since the random variables V and rj are un- 
corrected. One finds / [u,x] = 0, and 

f[u,x]f [0,0] ~ -A v d 2 x S(x) exp [-2ir 2 A a p \u\] . (7) 

Thus correlations in the pinning force f[u,x] are short- 
ranged in both x and u. [There is also an w-independent 
term (d x V) which is however irrelevant.] The dynamics 
of the random chain (5) (with the d x u term neglected) 
now look very much like the zero-temperature driven dy- 
namics of a directed path (DP), described by the trans- 
verse coordinates u(x), in an effective two-dimensional 
random potential ^[u,^] = — J u du'f[u', x] [3]. The 
convective term in (5) can be a relevant perturba- 
tion [9]. However, close to the depinning transi- 
tion, it is irrelevant since v — > 0. There, Eq. (5) 
taken on the simple form p,Q 1 d t u = — f^, with % = 
J dx [f (d x u) 2 + W [u, x] - Fu] , W = and 



W [u, x] W[u>, x'] ~ A v A a S(x - x')S(u - u'). (8) 

This result suggests that the critical behavior of the ran- 
dom chain may belong to the DP universality class, if the 
complicated correlations which exist in higher moments 
of the effective random potential are irrelevant [23] . 

To check our expectations, we performed numerical 
study of the bead-spring system. We first investigate 
the static behavior of the random chain in equilibrium 
using the transfer-matrix technique [15,17]. The system 
is simplified by discretizing the position of the beads. In 
Fig. 2 we present numerical data obtained from simu- 
lations of the system with TV = 4096 beads. Disorder 
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averaging was performed on 2000 different realizations 
of V(x) and a n . The potential V was chosen from a 
Gaussian distribution with zero mean and unit variance. 
The spring lengths a n were uniformly distributed in the 
interval [5, 15], and spring constants were kept equal at 
k = 0.2. From the preceding analysis, we expect the dis- 
placement field u n =r n — Ylm a ™ t° behave as the trans- 
verse coordinate of a directed path in a 2-dimensional 
random potential, with the total energy of the path E 
given by Eq. (1). Our numerical results for the fluctua- 
tions in the net displacement u n and the total energy £?* 
of the optimal path clearly indicate scaling (see Fig. 2), 



with ui ~ n 4 / 3 and AE* 



~ n 



1/3 



at zero temperature. 



Similar results are obtained at finite temperatures [15]. 
These scaling laws are characteristic of the 2-dimensional 
DP universality class [17]. 




FIG. 2. Fluctuations in (a) displacement, and (b) energy of 
the optimal path (measured with respect to a reference path 
which starts and ends at u = [15]). 

We next study the zero-temperature response of the 
random chain to an external driving force by numerically 
integrating the equation of motion (2). (Here, r„'s are 
kept as continuous variables.) The random potential V 
used is consisted of a series of bumps and valleys similar 
to those used in Ref . [24] . Each has a range b with an am- 
plitude uniformly distributed in the interval [— Vq, Vq\. In 
the dynamics, we explicitly put in an additional restric- 
tion which eliminates the crossing of beads [15]. This 
speeds up the dynamics without affecting the asymp- 
totic scaling properties. We consider open chains with 
free boundary conditions. The simulations were run for 
different chain lengths with parameters k = 0.1, b = 1, 
and Vo = 1. Various values of /io (time step size), rang- 
ing from 0.5 to 0.01, were used; we always checked that 
twice smaller /xo does not change results significantly. 

With a strong enough applied force F, the system 
moves and reaches a steady state after some time. The 
overall motion is defined by the steady state center of 
mass velocity v(F) = (v n (t)), where (• • •) de- 

notes here a time average for a given realization of the 
disorder. Averaging over long time is particularly neces- 
sary in the vicinity of the depinning threshold, F — > F+, 
since the motion of the chain there becomes very "jerky" , 
exhibiting "avalanches" in which some sections of the 
chain move forward with much higher velocity than oth- 



ers. Simulations were done for system size N = 1024, 
first for constant spring lengths a n = 10, and then for 
random springs, with a n randomly chosen from the in- 
terval [5, 15]. Depending on the value of F, time averages 
were taken over the intervals of 10 B to 10 7 steps. Simu- 
lations were repeated several times for different V's and 
a n 's; the average results for the velocity-force character- 
istics are shown in Fig. 3. The best power law fit to 
v ~ (F - F C )P gives F c ~ 0.275, (3 ~ 0.41 for the uniform 
springs, and F c ~ 0.289, (3 ~ 0.25 for random springs. 

The data presented in Fig. 3 clearly shows the differ- 
ence between the case of uniform and random springs, as 
anticipated from the above analysis. In the case of uni- 
form springs, our estimate for exponent (3 is comparable 
with the result (3 = 0.45 ± 0.05 obtained by Myers and 
Sethna [25] in their simulations of an automaton model 
developed to study the depinning transition in CDW sys- 
tems. Our result is also consistent with the FRG result 
of Narayan and Fisher [2] where (4 — e)-expansion yields 
(3 = 1/2. On the other hand, simulations of the random 
chain model give an exponent (3 which is close to the ones 
obtained in recent numerical studies of the driven DP in 
2-dimensional random media, where (3 = 0.24 ± 0.1 [6] 
and (3 = 0.25 ± 0.03 [7] were reported. 
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FIG. 3. Velocity-force characteristics v ~ (5F) 13 , where 
SF = (F — F c )/F c is the reduced driving force. 

In addition to the velocity-force characteristics, we an- 
alyzed the scaling of spatial correlations. Equal-time cor- 
relations between two beads separated by I springs are 
given by: 



C 2 (£,N) = N~ 



(9) 



We expect finite-size scaling to take the form C(£,N) = 
N^^(£/N), where $ is a scaling function. In order to 
estimate we simulated different systems from N = 8, 
through 256, close to the threshold F c . The chains were 
released from random initial conditions until they are 
(barely) pinned. Disorder averages were performed over 
many independent realizations, from 5000 for smaller 
iV's to 100 for the largest N. The insets in Fig. 4 
show numerical data for C (£, N) as a function of TV for 
fixed l/N . The slopes of these curves yield the expo- 
nent Q. To find the asymptotic exponent value, we cal- 
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culated the effective exponent at finite N according to 
C(N) = log [C(t,N)/C(%, f )] /log(2), and then averag- 
ing over all available values of I (see Fig. 4). The satu- 
rated (large TV) regime yields ( « 1.46 for uniform chains, 
and ( ss 1.22 for random chains. 
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FIG. 4. Effective roughness exponent C(-Y). 

Our value of £ for uniform chain is in very good agree- 
ment with the FRG result [2] of ( = 3/2 which is ex- 
pected to be exact. The result for random chain is also 
very close to the one reported by Leschhorn and Tang [7] 
for the driven DP in random-bond disorder (( ~ 1.25). 

Our numerical results combined with analytic argu- 
ment have demonstrated clearly that the critical dynam- 
ics of a 1-dimensional random chain in contact with a 
1-dimensional random substrate belongs to the univer- 
sality class of a directed path in 2-dimensional random 
media. We can thus expect the random chain to possess 
all of the scaling properties known for the DP, including 
the power-law distribution of avalanches, which means in 
this case a wide range of stick-slip events [26] . Our results 
on a random chain can be readily generalized to type (ii) 
systems in higher dimensions. Relevance of the latter to 
a variety of other problems including gel electrophoresis 
and pattern recognition will be discussed elsewhere. It 
seems hopeful that the disorder-dominated dynamics and 
statics of a generic type (ii) system may be understood 
by making analogy to the appropriate type (i) system. 
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